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PARTI. ANALYSIS OF A MULTI-STAGE AXIAL COMPRESSOR 


WITH HIGH REACTICN BLADING 





ABSTRACT FOR s:ART I 


An analysis of a multi-stage axial comzpressor with high reace 
tion blading is carried out. The usual methods of linearizing the equa- 
tions for flow through the cornpressor fail for this type of blading. A 
numerical solution of the non-linear equations is worked through for 
off-design operation. Flow velocities through the first six stages are 
calculated for the off-design flow rates in order to insure that the flow 
does indeed come to a steady-state repeating pattern. Entering flow 
angies for both rotor and stator are calculated for this repeating flow 
condition. The off-design incidence angles for the stators at low flow 


rates indicate a possible deterioration of efficiency. 





SUMMARY FOR PART I 


The pressure rise per stage in an axial flow compressor is pro- 
portional to the square of the incidence velocity relative to the blades. 
Almost all axial compressors for gas turbines were designed with the 
whirl in the same direction as the rotor blade speed. This choice of 
design flow gave high efficiencies when the blades were designed prop- 
erly. Until the introduction of transonic blading, a high efficiency was 
obtained only when the Mach number relative to the blade was less than 
unity -- usually in the neighborhood of 0.75. Relative velocities of this 
magnitude did not require whirl opposite to the direction of rotor blade 
speed and, in fact, the matching with a turbine was easier with whirl in 
the same direction as the rotor speed. 

An axial compressor for a gas of low molecular weight, for in- 
stance helium, can be designed with blading similar to the axial compres- 
sor for air. However, the pressure rise per stage is proportional to 
the square of the Mach number oi the relative velocity. The rotor blade 
rotational speed cannot be increased because of strength limitations, and 
hence the square of the relative Mach number and therefore the pressure 
yatio are decreased by a factor approximately proportional to the ratio 
of molecular weights. Thus, the number of stages required for a com- 
pressor using helium would be approximately seven times that for air 
for the same overall pressure. Clearly, higher relative velocities are 
required to decrease the number of stages for a gas of lower molecular 
weight. An obvious method of doing this is to introduce the whirl in the 
opposite direction to rotor blade speed. 

A vortex flow pattern was chosen as the design condition: then 
the flow is truly two-dimensional; that is, there is no radial shift of the 


streamlines as the fluid proceeds downstream through the compressor. 





This considerably simplifies the analysis of the flow at the design flow 
rate. Off-design calculations must still be carried out in order to in- 
vestigate the compressor performance and to indicate how rapidly losses 
might increase as the compressor flow rate increases or decreases from 
the design value. For conventional multi-stage compressors with rela- 
tively low ideal pressure coefficients, this analysis is usually carried 
out by linearizing the governing equations for three-dimensional flow 
through the compressor, and is quite satisfactory. However, with whirl 
against the rotor blades and its correspondingly higher pressure coef- 
ficient, this type of linearization fails completely. It is necessary to al- 
low for the streamline shift to avoid singularities in the solution. Thus, 
the full set of non-linear equatinns, taking into account the radial shift 

of the streamlines, must be used. In order to calculate the repeating 
flow pattern characteristic of multi-stage compressors, it was necessary 
to numerically determine the flow pattern through the first four or six 
stages. 

Although detailed calculations of efficiency are scarcely possible, 
the results of the analysis can be used to give indications of probable 
losses. The degree of turning required by the blades appears to be with- 
in permissible limits. The off-design incidence angles for the rotors 
deviate very little from the design values. The stator incidence angles 
are more severe, and might produce a loss in efficiency over part of the 


stator blades. 





I. INTRODUCTION 


Conventional subsonic axial-flow compressors designed to use 
air as a working fluid are not necessarily sitable for gases of lower 
molecular weight. For a particular closed cycle system, where a 
working fluid other than air might be desirable, the number of stages 
required will vary greatly with the choice of gas. The number of stages 
is directly proportional to the specific heat at constant pressure and is 
inversely proportional to the blade speed and the change in tangential 
velocity of the gas through a stage. The blade speéd is fixed by the 
strength of the blade material and is generally held as high as possible. 
With present-day air compressors, there is a maximum amount of 
tangential velocity change that can be accomplished efficiently. Thus, 
for a typical compressor designed to use air, the number of stages re- 
quired would be increased approximately five times if a gas such as 
helium were used and the same temperature ratio across the compres- 
sor were to be maintained. This report describes a blade design that 
will hopefully reduce this high number of stages required for helium to 


approximately one-half the number presently used. 





Unlike conventional multi-stage compressors, the design of the 
stator blade is such as to give an accelerating flow through the stator 
rows. The rotor blades then turn the flow back to the axial direction. 
This type of design has been used in some single-stage fans. The nota- 


tion to be used in this report is illustrated in Figure 1. 
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Figure 1, 


The following list of symbols will be used throushout this report: 


c.. * axial component of velocity 
cy * tangential component of velocity 
c¢ = absolute velocity 


w velocity relative to the rotor 





o3@ 
=, = tangential velocity relative to the rotor 
y 2 angle between axial direction and absolute velocities 
3 = angle between axial direction and relative velocities 
R = tip radius 
us =-—stip speed 
Ww = angular velocity of rotor 
r = radius measured from compressor centerline 


static pressure 


*S 
iW 


9 = density (assumed constant throughout) 
Subscripts refer to stations as shown in Figure 1. 

In two-dimensional" flow through an axial compressor, the 
streamlines are assumed to lie on cylindrical surfaces with constant 
radii. The pressure gradients in the radial direction are ignored, al- 
though it is shown later that at design for vortex blading, the flow is 
two-dimensional and the pressure cradients are properly accounted for. 
Since an important difficulty in applying two-dimensional flow theory to 
this design arises, we will discuss it briefly here. 

The blades are of the free vortex design. Thus 

“y, “« l/r. 
A dimensionless radius is defined by 
= r/R 
Therefore, the tangential component of velocity through the stator is 
taken as 
“a, = ua/E , (1) 
where a ic an arbitrary constant to be determined. 


The pressure drop across a stator row is given by 





a. = 


} So £ at Z 2 
eee Blinc? -c,*) . (2) 


Since r is a constant along a streamline, the Bernoulli cquation for 


rotating axes sives, for the pressure rise across a rotor row, 


i aT } 2 c 
a a = AW, ae ) e (3) 


The pressurerise across a stage is then 
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At design, oc. = QO, and from the velocity diaxvrarns we obtain 
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a 
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After substitution of these relations into (4), we obtain, for the pres- 


sure rise across the stage, 


= uo , (5) 


A work coefficient is usually defined as 


y = = . (6) 


spw R 
Thus, when the compressor is operating at the design point, 
" = an. 


However, & is not unrestricted. It has been found experimentally 


desirable (Ref. 1) to have 





——s-~— <5. (7) 


where 8 is between 9.5 and 0.6. Defining a non-dimensional Mow 
coefficient as 


¢p = ci /u : (8) 


this restriction can be written 
2 
e526) : .s lee (9) 
af ~ {1-8 &~/6 

This condition must be satisfied over the entire length of the blade. 
Since the restriction is most severe at the hub, values of § and 4 
evaluated there will be conservative over the rest of the blade. An ex- 
isting experimental set-up prescribed the hub ratio Fr =90.6. There is 
another restriction on the velocities that is not easily degcribedinanaly- 
ticalform. This is the amount of turning that can be accomplished ef- 
ficiently in the blade rows. Due to lower blade speed at the hud and the 
adverse pressure gradient in the rotor, this condition also is most 
severe at the rotor hub. 6 was chosen as 0.6, and values of the turning 
angle were obtained for different values of ". The quantities ' = 0.70 
and ¢ = 0.53 were subsequently chosen. The value of ¢ can be raised 
substantially for a compressor with hub ratio in the range 0.7 to 9.75. 

So far, no difficulty with the two-dimensional theory has arisen. 
This is not unexpected since, as montioned before, vortex blading is 
exactly two dimensional at the design flow rate. However, we have yet 
to consider the off-design conditions where a radial component of ve- 
locity does exist. At design, the work coefficient was constant over the 


blade height. This is desirable in order to cut down losses cue to mix- 





-6« 
ing at the end of the compressor, and actually is neccssary on all stages 
after the first few because a piling up of pressure at sore radius is 
impossible for a very large number of stages. Thus, this restriction 
will be applica at off-design flow rates also. Then 
vy = ad = mT (w “aw, 405 “e0,°) = constant. (10) 
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r=) 


J 


Remembering that c # 0 for off-design conditions, then 
fa 


w = ub +c : 
a bee | 


ae = uf + en. ry 
z to 
and the work coefficient becomes 


& 
y= 22 (c -c ) 
u | Us 


w = 2ele + o(tany, - tan3,)) : (11) 
A. first approximation can be made that tany, and tan B » correspond- 
ing to the leaving angles of flow, are independent of the inlet angles and 
hence are independent of flow rate. If this is done, the local value of ¢ 
at the root is 


' : Ori? = 


Since under our approximation { tans, - tany , ) is constant and is great- 
er than zero, we see that for * > 0.72 the flow coefficient at the root 
becomes negative and is meaningless. If tan’ , and tany, are allowed 
to change with flow speed and this change calculated by cascade theory, 
the singularity in is not relieved appreciably. Since values of the 


ideal pressure cocfficient greater than 0.72 would be required in starting, 





ay 
for instance, the two-dimensional analysis of the flow is unsatisfactory. 
It must be mentioned that if + could be pushed as high as 0.72 ina 
standard compressor, there would be the same difficulty. Normally, 
axial-flow air compressors have an ideal pressure coefficient of be- 
tween 0.4 and 0.5. It will be shown that a three-dimensional flow 
analysis of the compressor does not contain this singularity in 6. 
This difficulty is mathematical rather than physical, and appears to be 
the result of an invalid linearization of the three-dimensional equations 
of motion. There seems to be no way of calculating the ideal flow pate 
terns except by a numerical procedure with a radial shift of the stream- 
lines taken into account. 

The assumptions made in the three-dimensional flow analysis 
are: 

(1) incompressible perfect fluid; 

(2) axisymmetric flow; 

(3) equations are written for conditions far upstream and far 
downstream of a blade row. 

The notation used is the same as in the two-dimensional analysis 
except that the odd subscripts refer to stations far downstream of stator 
rows and the even subscripts refer to stations far downstream of rotor 
rows (see Figure 2). 


The Bernoulli equation for flow through a stator row is 
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The Bernoulli equation for flow through a rotor row is 
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Figure 2. 
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The equations for the pressure gradient in the radial direction are 
. 2 
| Cc g 
Pan | 5 2h ang Pamth 2 | Mane (5) 
= *2n ae haa *2n+1 


By using (15) to eliminate the pressure terms, equations (13) and (14) 


become, respectively, 
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Considering a cylindrical stream tube, the continuity equation can be 
written in the differential form 


dr : 


Cc. dr = <¢c s 
Ao n+] 2n¢1~° 2n+i 


a, 2 on 


Within the scope of the assumptions previously mentioned, these equa- 
tions are quite general. It is now assumed that the leaving angles of 
both stator and rotor are constant with respect to flow rate and give 
vortex flow at the design flow rate. Thus, evaluating the leaving angles 
at the design condition, indicated by subscript d, we obtain 


C, 





Cc =z C tany = 9 (18) 
Vent] arnt] ee 
d d 
or, in the non-dimensional notation, 
tany, = = oe : {19) 
1 Oa ont) %9a Song) 
Similarly, we obtain for the rotor leaving angle 
S 
imep, SSD . (20) 
94 


With the leaving angles so specified and the following definitions, 


Pl”, ie. and b = 9a. 


the equation for the flow through the stator becomes 
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The equation for flow through a rotor row becomes 


do on 
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The continuity relations are 


anelSan-1"S ane = P2nSan°ban = P2ntiS2nti Santi > (23) 


The two pairs of non-linear equations (21) and (22), and (23) can be 
solved approximately for 6, 38 a function of Em fom mt 0, 1772,.2. 
in succession. 

Notice the particularly simple solution by? Oy for all values 
of m. It follows that there is no streamline shift and no radial velocity 
in this case, and the analysis of the flow reduces to the two-dimensional 
one previously considered. It is largely because of this analytical sim- 
plification that the free vortex design has been used so frequently. How- 
ever, for off-design flow rates, no such simple solution exists, even 
for the free vortex design. 

To illustrate the method of solution to these equations, let us 
consider an off-design constant axial flow 6 approaching the first sta- 
tor. In this case, care must be taken in applying equation (21), since 
the assumption of a whirl velocity is incorporated in the right-hand side 


of this relation. [For this incoming flow, equation (13) becomes 


‘ 


Zp 2p > 
p l l P 0 


Since dp /dr, = 0, equation (16) becomes 





Z2 
a, , 
j d _ 
= + aS (c., “- Se, =e 





alee 
Substituting from equation (19), the velocity diagram relations, and the 
definition of a, the equation corresponding to equation (21) is 
2 d u 
(1 + : x 0 
pel dé 
ilies 
or 
>) = constant. 
This constant is then determined by integrating the continuity relation 
from the hub to the tip. 
. l 
J $696) = | ae 
S A 
Z é 
oy z(1 - 6.) > % ail = &, ) 
or 
) a Po : 
Equations (21) and (22) are applicable to all of the other blade 
rows now. The relation for > becomes, using equation (22), 


5 Z 
J b +6, 
Ce rr Soe 


b 
Integrating, we obtain 
nae? 
2 


1 ae oe 
D +§&5 


where K is determined as before by integrating the continuity expres- 
sion. 

The flow coefficients for the following blade rows are obtained in 
the same manner; however, the integration for these must be done nu- 
merically. Also, another mathematical difficulty arises in the solution 


for $, and the remaining stages. The existence of both independent 





variables on the right-hand side of both {21) and (22) will reguire an 
iteration process. The procedure is to first let b ont] = 5 on and solve 
for don 41° The continuity relation can now be integrated from om to 
€, and a relation between bon 7" and € 2n obtained. Substitution of 
this relation into the differential equation for Pong] yields a refined 


value for 6 Fortunately, the correction to the flow coefficient is 


2n+1 ° 


not large, and the iteration process is senerally not necessary. 
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Il. APPLICATION OF THE THREE-DIMENSIONAL THEORY 


The calculation of the flow field as it progresses down the com- 
pressor was carried out for two off-design conditions. The blades were 
chosen such that by = 0.57751 and oe = 0.70000. The average flow 
rate for the off-design conditions was taken as $ = .46000 and ¢ = 
.69000. Figures 3 through 6 show the variation of the flow coefficient 
versus dimensionless radius at various stations in the compressor. 

The flow field at lower than design flow rate is converging fairly rapidly 
to a steady-state condition, while at the higher flow rate the convergence 
does not appear to be as good. 

The error in choosing bond] = bon as a first approximation is 
shown in Figure 7. The calculation was carried out for the fourth rotor 
row. An error of only approximately 2 percent of bg at the mid-height 
of the blade is introduced by the first approximation. It is felt that real 
fluid effects will be greater than this amount. 

The local values of " for the two off-design flow rates are shown 
in Figure 8. Notice that ! is very nearly independent of the radius, 
agreeing with our assumptions leading to equation (10). 

In standard compressors there esadsts a perturbation solution to 
the three-dimensional equations. There it is assumed that the flow om 
from a stator ic approximately equal to the flow 9 from arotor. Thus 
the procedure is to find a solution for which loo, | << ($5+0,). To 
the zeroth order, this results in a linear compressor performance 
curve in the form 


h# = A-BO. (24) 





-]4- 
The two-dimensional theory also produced an equation of this type, 
since in that case $y was taken identically equal to 5 . 
To examine the possibilities of this approach for the proposed 
helium compressor, consider the approximate performance equation 


? 
«* = 2fE~ + Ed tany, - Eo,tans, | 2 (25) 


For free vortex blading, equation (25) becomes 


te 2 ae ere ; (26) 
If we now assume 
b, = 216) +62) + 21,762) = O, + Ad (27a) 
and 
$5 = 2, +5) - 2(9,-$,) F 6, - AO, (27b) 


then equation (26) can be written 
ae 24. Ad 2 
# = 2€ ae ag G+ ee ltgt ee) - (28) 


Typical values of #, and Ad for a conventional compressor allow the 


a 
term containing Ad to be ignored. However, for the compressor pro- 
posed in this report, Ws is approaching ae ° at the hub, and the term 
containing Ad can no longer be neglected. Figure 9 shows the per- 
formance curves of both types of compressors with and without the in- 
clusion of the Ad term. The inadequacy of the two-dimensional 
solution in the case of the helium compressor is obvious. To the scale 
of the graph, no difference can be observed for the conventional com- 
pressor. 


Figure 10 shows the total pressure variation along the blades 


after the flow has essentially settled into a steady pattern. For lower 





Sie 
than design flow rates, more work is being done by the tips, while at 
higher flow rates, the roots give the greater rise. The cumulative dif- 
ference after six stages is still less than half the design pressure rise 
for a single stage. 


It was mentioned previously that it is desirable to have 


ae ae.” 
\ é < 0.6 


Wi l 


This quantity is plotted in Figure 11. FPossible difficulties are indicated 
at lower than design flow rates. A higher hub ratio would eliminate this 
difficulty. The corresponding ratio for the stator is plotted in Figure 
12. Since the stators are in a field of decreasing pressure, no diffi- 
culties are expected here. 

The off-design entering angles for the stator and rotor are plot- 
ted in Figures 13 and 14, respectively. The stators are very critical 
and will require care in their design. Sharp leading edges of the sta- 
tors would probably cause high losses. The off-design flow incidence 


angle deviations for the rotors appear satisfactory. 
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IV. CONCLUDING REMARKS 


It has been demonstrated that high reaction blading can be used 
in multi-stage axial compressors to produce a sipnificant increase in 
pressure rise per stage over conventional blading. This type of blading 
could be used in a closed-cycle helium compressor and greatly reduce 
the total number of stages required. ‘The analysis also indicates that 
a hub ratio greater than 0.6 would possibly be a more efficient design, 
since stator incidence flow angle deviations become quite large near the 
hub in the example shown. 

It was also shown that the usual linearization of the three- 
dimensional equations of flow is not valid for this type of blading. How- 
ever, a numerical solution allowing for the shift of the streamlines does 
converge to a steady, repeating flow pattern after the first five or six 


stages for flows as much as twenty per cent off design. 
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Figure 3. Stator Axial Velocity Coefficients, 6 = . 46000. 
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Figure 4. Rotor Axial Velocity Coefficients, $6 = .46000. 
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Figure 5. Stator Axial Velocity Coefficients, $ = .69000. 





ST hes 








QO400 
0.6 O07 0.8 Oe) lO 


: 


Figure 6. Rotor Axial Velocity Coefficients, $6 = .69000. 
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Figure 7. Streamline Radii Upstream and Downstream of a Rotor. 
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Figure 8, Local Value of ¢ Versus Non-Dimensional Radius. 
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Figure 10. Total Pressure Variation with Radius. 
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Figure 12. Kinetic Energy Conversion Through the Stator Versus 
Radius. 
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Figure 13. Stator Flow Entering Angles. 
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Figure 14, Rotor Flow Engering Angles. 
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Figure 14 (continued). Rotor Flow Entering Angles. 





PART II. A DESIGN STUDY OF A MULTI-STAGE AXIAL 
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ABSTRACT FOR FART I 


A method for designing blades of high aspect ratio in a multi- 
stage axial compressor is investigated. This method requires a cal- 
culation of the flow induced by all blade rows. The method is illustrated 
for the limiting condition of hub ratio equal to unity and then is carried 
out for the more realistic case of finite hub ratio. 

An example of a blade design for a particular flow is carried 
out. The results are compared with those of a previous design for the 
same flow conditions but based on a theory applicable for blading of low 


aspect ratio. 
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SUMMARY FOR PART Il 


The flow through most present-day axial flow compressors is 
frequently analyzed by applying a three-dimensional theory that is based 
on the assumption of blading with low aspect ratio. Even though existing 
compressors do not have what could be called low aspect ratio blades, 
the theory predictions and experimental results show reasonable agree- 
ment. 

Where weight considerations are of primary concern, the possi- 
bility of reducing the size of a compressor by shortening its length is 
attractive. There is no known reason why the blade rows cannot be held 
to their same height but shortened in chord. There is, however, a dif- 
ficulty in applying the aforementioned analysis. There, it is assumed 
that equilibrium flow patterns exist at the leading and trailing edges of 
the blades. This is identical to assuming low aspect ratio for the blades, 
since they must be wide enough to allow all changes in the flow to occur 
entirely within the blade row itself. Several extensions of the theory 
have been developed to take account of the non-equilibrium flow near the 
blade rows. Some of these analyses assumed small vorticity throughout 
the flow field and also a small change in vorticity through the blade rows. 
Another theory relaxed the restriction of small vorticity throughout the 
flow field. 

The theory presented in this thesis is an extension of the latter 
approach, and is based on the assumptions of closely spaced blades of 
high aspect ratio. The equations of motion are linearized by a perturba- 
tion scheme that assumes small changes in perturbation velocities across 


the blade rows. To illustrate the method, the calculatinn is first 





~ae 
restricted to the limiting case of hub ratio close to unity. Then the 
more involved problem for any hub ratio.is investigated. An example 
is worked for blades of aspect ratio about three. Here, there appeared 
little difference between the theory applicable to blades of low aspect 


ratio and the theory of blades of low aspect ratio. 
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I. INTRODUCTION 


The theory of three-dimensional flow in the axial compressor 
described in Part 1 of this thesis is concerned only with flow conditions 
far upstream and far downstream of a blade row. This theory would be 
applicable to blades of very low aspect ratio. The difference between 
these velocity profiles and the actual profiles near the leading and trail- 
ing edges of the blades was assumed to be negiigible. It was on this 
basis that the entering flow angles for the blades were determined. If 
the aspect ratio is not very low, it is possible to obtain more exact in- 
formation about the flow pattern near the blade rows by making use of 
some of Marble's results (Ref. 1). This has been done by Bowen, 
Sabersky, and Rannie in Ref. 2. The method carried out by these in- 
vestigators forms the basis of the theory to be discussed in this report. 
The approach is to consider the disturbances due to rows of blades of 
high aspect ratio placed very close together. Here, the flow in any par- 
ticular blade row is influenced by neighboring blade rows. In the three- 
dimensional theory of Part], the changes in the flow pattern were as- 
sumed to occur entirely within the blade row itself and one blade row did 
not affect its neighbor. 

The analysis begins with the calculation of the disturbances due 
to an infinitesimally thin disk. These disturbances can then be super- 
imposed to form a blade with finite axial extent. The effect of many 
blade rows can be determined by the superposition of the disturbances 
due to the finite row. In this way, the flow pattern for a multi-stage 


axial turbomachine can be constructed. For a repeating flow pattern, 


the summation of the effects of all blade rows is possible. 
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i. TV'O-DIMENSIONAL THEORY 


A. Equations and Boundary Conditions 


Before proceeding directly into the three-dimensional theory of 
high aspect ratio blading, it is informative to investigate first the limit- 
ing case of high hub ratio. This problem is simpler because the radial 
pressure gradients vanish far from the blade. The analysis is based 
upon the disturbances due to a disk at x = 0 (see Figure 1), and the 


equations of motion are written subject to several assumed conditions. 


pr 





Figure | 


These assumptions are: (1) there are no circumferential variations; 
that is, = ()=0. This is equivalent to having a disk with an infinite 


number of blades. (2) The flow is incompressible, the density p being 
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held constant. The continuity and momentum equations then become 


mtg: 0 (1) 


Ou, Su _ .} 
” Ox ay 2 
ee. y wy 1 oO 
tx + % By 7 5 


eho 


(2) 


(3) 


we 


9) (4) 


3 
a 

“+ 

< 
Si 
ii 


where u,v, and w are the gas velocity components in the x, y, and 
Z Girections, and p is the static pressure. 
The boundary conditions to be imposed on this system are as 
follows: 
(1) vz0 at y=0 and y=b forall x, since the walls of the 
compressor are streamlines. 
(2) u is continuous at x= 0 from continuity through the actu- 
ating disk. 
(3) v is continuous at x = 0 due to the component of force 
being zerointhe y direction. 
(4) w is allowed to change discontinuously across x=0@. 
(5) At x=s-o, 
u = uly) - cugty), 
v= O, 
w = Wy) = ew ly), 
p=p- ep, = constant. 


(6) At x=+ 0, 


iF 


a uly) 3 cu ily) 9 


vee 0 , 
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w = wy) + ew ly). 


p = Pp + Ea = constant. 


These last two boundary conditions are written in a form so as to con- 
form with the following perturbation scheme. Far upstream and far 
downstream of the disk, the flow is independent of x; any disturbance 
due to the disk having been smoothed out. Close to the disk, the flow is 
disturbed by a small amount, and the following perturbations are as- 
sumed, where the symbol e is carried only to denote the assumed 
small quantities. 

ui = uly) + cuply) + cu, (sy) 

v = 6 V(x y) 

Apeicasin (5) 

w = wy) + ew {y) + ew (x y) 

p = pep, + ep, (xy) 
Notice that p and Py are assumed independent of either x or y and 
that the mean value of v is zero. The top signs are for -am<x<0 
and the bottom signs are for 0<x<+o. 


Introducing equation (5) into the equation of motion yields the 


linear set 

Ou, avy 

ax tap 7 (6) 

ae oe (7) 
Ox oy i: o ox 

- ov, nee oP, ie 
Ox p oy 
Ow = 

= l Ow a 

ae ey iy (9) 
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B. Solution of the Equations for the Actuating Disk 


The first three of the preceding equations (6, 7, and 8) can be 


combined to give the differential equation for v, as 
2 


2 
Rig ta ey an 
<r e. ee ae ey 
Ox oy u dy 


with the boundary conditions that v=0 at y=0 and y=b. A solu- 


tion to equation (10) that satisfies these boundary conditions is 


+i : 
v, = U y AL e- By Wy) (11) 
where 
b 
o¢e laa 
B. y 
0 


and YAY) satisfies the differential equation 


(@)-2]s, 


with boundary conditions 


a 


ti 


0 


Y,(o) = ¥_(b) = 0. 


Equation (6) then yields 


sy 


= +d 
a4 © = ta B viyyy, (12) 


where the prime denotes differentiation with respect to y. Notice that 
the boundary conditions on u are satisfied at x=-o and x=+o. 
rom the continuity of u across the disk, the following relation is 
obtained 


CO 
uly) = - U ) AL = ¥'(y) « (13) 
n 


n=l 
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Since u and u, are to be chosen and the shape of the blades is to be 
determined, it is possible at this point to choose AL =OQO for n>l. 


Also, if u is chosen as a linear function of y, then 


¥y = sin a 
and 
hy = 7 . 
Note that it is not necessary that these particular choices be made. It 
is done here only to reduce the amount of algebra in the rest of the analy- 
Sis. It will be pointed out later that these choices can indeed represent 
the flow through a conventional compressor quite well. Therefore, the 


perturbation velocities become 


+uE 
a, od ay nad 
and 
tug f 
¥, = -@ | u(n)dn . 
Ww 0 
0 


Using equation (9), Ww) is determined as 


x 
— +7 
4 ] di \—F5 
w sto e i dn, (14) 
> 0 
and with the use of equation (7), 
tmE da F 
l _7- aad u 7 
P rq = + €e E os “dy , wi(n) én. (15) 
0 


There are several requirements that we can now impose. First, 
we can require the total pressure rise across the disk to be a constant. 
This will insure that the flow pattern will repeat itself. This is a similar 


condition to that set upon the work coefficient in Part I of this paper. To 
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calculate the change in total pressure, one needs the change in whirl 
and the change in static pressure. The change in whiri is 

Aw = w(to,y) - wl-o,y), 


and with the use of equations (5) and (14) becomes 


/ 
9 


| 


<3 


7 
aay 
E. 


V 


W = 2cW. = 2e i. 
rAN 7 u(r) dn 


The change in static pressure is 


Ap -P FE X.. 
r ‘ (+0, y) : (-0, y) 


and using equations (5) and (15), becomes 
a : 24 vlsu fy ine 
pa fo 7 els @ Oy Jeo | * 
0 


Since u and v are continuous through the disk, the change in total 


pressure is 
Sp 
ee q = = +w Aw 


or 


AP. Pp i « : x 


and we require this to be constant. Thus, we now have the four unknown 
functions of y: u, use w, and Woe and two constants p and Po 
with one equation (16) relating them. Also, since we are interested only 
in the pressure rise in a compressor, we are not immensely concerned 
with the mean pressure, p. The pressure change across the disk is 
25 » and this is the more important quantity. 

Another condition that we can require is that the total pressure 


far upstream and far downstream be constant from root to tip. Note 


42-6 
that this is not a necessary condition, but is a common choice in come 
pressor design and gives two more relations connecting our four un- 
known functions of y. For the total pressure to be constant we must 


require that 


constant 


Oe Pinu.) 4 (arten.)- 
= + sf{(u eu.) + (w ew.) 1 
and (17) 


5 +3 [(Bteu,)* 2) (Few) J = constant 


and since p/p is constant far upstream and far downstream of the 
disk, we derive the following two equations by adding and subtracting 


equations (17) and keeping only first order in « : 


a° sa" = Ky 
_ ~ (18) 
Got wy = Ks 


where Ky and Ko are undetermined constants. 


The second condition of equations {13) is now actually the same 


as equation (16). If the disk were to act as a stator, then (16) would 


reduce to 
22% +W_We- GTA) fy ye (19) 
eee: et (3 = @/ Je 
0 


Since p,/p is a constant, and from (18) uuij t+wiwsK,, the last 
term of equation (19) must be a constant. Since the integral is obviously 
zero at y= 90, then this constant must be zero. Thus, with the assump- 
tion of constant total pressure, we have four unknown functions of y 

and two conditions upon them. Hence, for example, we could pick u, 


use and Po? and the flow will be completely determined. Once this 
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ew béen cCeme, the shape of the blades can be found. 

To derive the effect on the flow cue to a blade of finite width a 
and height b, we can sumfrom x= 9 to xza the disturbance due to 
a disk located at x =§ where § rangesfrom 0 to a. Notice, how- 
ever, that here one must make a choice of the shape of the loading on 
the blade. That is, we must require that each disk contribute some 
fraction £(&)c& of the disturbance of the blade. In the problem pro- 
posed in this report, we are at liberty to pick the function f£(€). We 
then go on to find the shape of the blade. However, the inverse prohlem 
of specifying the blading and then finding the flow is much more diffi- 


cult. In that case, the loading function is in general a function of both 


x and y and depends on the entering flow angles. 


C. Superposition of Actuating Disks for Multi-Stage Compressors 
Given the function f(£), we can proceed to determine the ve- 
locities u and w throughout the compressor. The function {(€) is 


normalized such that 


For a finite stator of width 2a, height b, and leading edge at x=0, 


the axial velocity for x>0 becomes 


x a x 
&. = 7 “a ee 
us ue cu, + eu, e€ | 6) e de. (20) 
rer 0O< & < #, 
P< 6 a Tt 
a a - ~(x-§) + w= (x-6) 
uszuteu, "ag)|1- 6 b Jat - aig) [1 =e o Jae » (21) 
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For x>a, 


-Tr 7 +1 
us uteu- eu, e | 46) e dé. (22) 
0 


The equations for a rotor are similar with only the algebraic sign of 
Uo changed throughout. 

Equations (20) - (22) give the flow due to one blade row alone. 
For a repeating flow pattern, the velocities are the same as if there 
were an infinite number of blade rows. Thus, for the complete flow pat- 
tern, we must superimpose the flows due to many stators and many ro- 
tors. The following definitions will be used: 


width of a stator row 


rs} 
it 


width of a rotor row 


@ 
FT 


Q 
WW 


gap between stator and rotor 


gap between rotor and stator 


ie) 
ti 


£ (6 )a—E = fraction of disturbance due to stator 
£6 )d§ = fraction of disturbance due to rotor 
4 = ato, Pate. 
The leading edge of the stators are then set at 
x = mt where m=0, +1, +2, ... 
and the leading edgesof the rotors are at 
x = a, te, tmet , 


The axial velocity in a stator row due to all blades becomes 
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x wT ar ct) 
a =p e-8) +H (x6) 
u = Uly) + cu, (y) feel -e» lee . Pee -e » Jal 
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0 
a 
9 6) r T 
ra) 'Xe-a =-C -mt ] 
re - unfy) | £,(6) e wt, gf ¥ 
m=2Q 0 
a. 
- ange [xem ] 
+e a u(y) See « € P e B 
m=] 
= “e -1 S + [x-mt] 
+e » u(y) | £,(6) e 5 ap ec 5 
m=] 0 
a 
= - tae =o [x+m2] 
-e uty) | Be ab | (23) 
ml 0 


Carrying out the summations, equation (23) becomes 
x T * 1 
ie - 5 (x8) ; e + 5 (¥-S) 
a= By) + eugty) J #6 l-e ae - | £8)|1 - 6 ag 


€ on > T -~(a_tc_) 
* | foo Ea 2 Seb mies 


lee 


E x 
=. 5 (a,+¢,) ~T Te oe t 
we” a BSB fee Page e © 


a 
S 
r tre “TF “Bef 


+ £ (8) e dé e e (24) 


oc 


For the axial velocity in a rotor row, change the sign of u * and inter- 


change the s and r subscripts in all terms of equation (24). 





The corresponding calculation of the whirl velocity in a stator 


results in 


Pp j 
o z = 5 (x-§) 
wetyte [£G)lwoy)-+ 2 | ainane ® |e 
0 . 0 


a 
S 1 
p = + x (x-§) 
an f (5) w (7) - zs — uolr)dn e 5 |e 
x i 0 
x Oy E = 2 
es -7 a -y-(a_+c_) 
ci = i. a (n}dn ———ay = fe, (5S) e Pag e Fy FS 
: 0 \ac. 
-¢ 
at € - . ae x 4 
+1 -7 e~(a_tc_) <7 
fatere a a eee "E een 
0 
a 
‘ + 7 swe Sac 
+ le (&) e fas o, 8 : (25) 


0 
The whirl velocity in a rotor is determined by changing the sign of a. 
and ss and interchanging the s and r subscripts in equation (25). 


An example will be shown later for which the following assump- 


tions are made: 


£(é) = £(€) = l/a, 


L = 2#(ate). 


D. Flow Patterns for Blading of High Aspect Ratio 


If then the exponentials are expanded in polynomials of a/b, 


c/b, or x/b, as the case may be, and terms through second order are 





a, 2 


retained, the velocities reduce to 


us Uteuw(F-) 
(26) 
we WHew (25-1)- fF (22-1- 58458 ) 
for a stator, and 
us U-runr(s- a 
a = . xo x* _ Tx Ta 
w= ew, (25 L) tet (er ] = ye? 
for a rotor. In both cases 
1am f 
[ £ ay | ui(n) dn : (28 ) 
? 0 


The entering and leaving flow angles for both stator and rotor 
blades can now be calculated. The defining equations for these angles 


are,in the case of a stator, 


_ u(x = 0) 
fan ¥) = wie s Oy 
and (29) 
_ u(x =a) 
ae. - wie & 


where Y) is the etator entering angles and Y> is the stator leaving 


angle. For the rotor, 


u(x = 0) 


tear) = Glade speed - w(x=0) 


u(x = a) 


tan P> = Blade speed - wixea) 
In summary, given the blade loading, blade aspect ratio, the 
axial velocity through the compressor, the whirl velocity at one radius 
of a blade row, the desired pressure rise across a stage, and the de- 


sired degree of reaction, one can compute the shape of the blades subject 
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WI. THEORY FOR FINITE HUB RADIUS 
The analysis for the limiting case of hub ratio equal to one is 
extended to the axisymmetric case of finite hub radius with only slight 
modifications. The y coordinate goes over into the radius r. The 
velocity inthe r direction remains v. The major change here is that 
the radial pressure gradients must be accounted for properly. 


The equations of motion in the axisymmetric case are 


Ou, 
dr * = i 
ou * 1 ap 
be a ieee p ox 
(31) 
. Ago ov  w _ _1 Op 
- Sr > * "ptr 
Ow ow “ 
u Bo + Vv ’ Fr ti £6 0 
As before, we assume a perturbation approach of the form 
us u(r) + eu (r) a EU, (x r) 
w= rv (X r) 
(32) 


w = wr) + ew (x) + OW y (Xs r) 
P = P(r) + ep i(r) + ep, (x Fr) 
Notice that p and P, are now functions of the radial variable, r. 


The resulting linear equations are 
ou 


1,19 a 
ox tr Or (FY) = 
= oud PY 

“Ox 10r — 0 Ox 
+ A 4, 8Py (33) 
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and also the radial pressure equilibrium equations 





om ie, 
1 oO WwW 
——— > (34) 
and _ 
1 Sp , ; ew WO 
op Or ~ r 


which must be satisfied at x = +o. 


Eliminating the other dependent variables, the equation for ay 





becomes 
av, Ov, 1 oy 
a t= e - (st Plt) lv, = O, (35) 
x r 
where 
= rena. ,h God 2w d 
Bi) sae (fF ae) - 9 te oe oe) 
u ru 
An appropriate solution is of the form 
HAE 
Vv, = * R(r) . 
Substituting into equation (35) yields 
ri 42 ate: : F(r)jJR = 0 (36) 
Fa Ue? taco ee 


r 
where the prime indicates differentiation with respect to r. If F(r) is 
chosen identically equal to zero, then (36) reduces to Bessel's equation. 
If we then take only the first term of the solution for Vio the boundary 
condition that v= 0 at the hub and the tip gives, for the first eigen- 
value 

hy = 3.18 
for a hub ratio of 0.6. If F(x) is small compared to (x, /e)° » then 
the first cigenvalue does not change substantially. We will now choose 


F(r) identically equal to zero and later compute its actual value to 





«ies 
: 2 
compare with (,,/b) : 
Proceeding as in the case with infinite hub radius, the perturba- 


tion quantities become 


the (37) 


whe re r, is the hub radius. 
The condition that the total pressure rise across a disk be a 


constant yields 


r 


eo O _ lduw ,dw.w C 
ee = et 2-(au,+Ww,) - 22+ Z (E+E) | NBy(nidn.(38) 

r 

i 


Again choosing the special condition of constant total pressure in the 
radial direction implies to first order in « that 
Dt+3p (7 +u’) = Ky 
and (39) 
Pote (ww + uu) = K,- 
Then, as in the simpler case, equation (38) applied across a stator re= 
quires K,=0. Combining equations (34) and (39) results in the two 


conditions 





ied eo. oe 
= 5 Amarr 


and (40) 


= (2° Fw) ors = (uu) , 

The remainder of the analysis, consisting of the summation of 
the disturbances due tc the disks in order to produce a finite blade and 
the summation over all blade rows to derive the complete flow pattern 
for the compressor, proceeds exactly as in the limiting case, and only 


results will be given here. ‘With the definition 
f 
dw 
r = (a+ Z) [mucin an. 
*t 
the velocity components for the stators are 
are x a 
us ut eur (= a 
and (41) 
S ue ae x. - X_, Ax 4 he 
ws w+ ew (2 7 1) cP (2 "i ] = +35) 
and for the rotors 
ie x a 
u= Ue eur (= -sr ) 
and (42) 


The blade angles are determined exactly as before with the use of equa- 


tions (28) and (29). 
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iV. AN APPLICATION OF THE THEORY 


Before proceeding with an example that will allow us to verify 
this theory against an existing compressor, we will indicate the steps 
in order. 

1. Choose the blade loading function {(£) for both the rotor and 
the stator. 

2 Choose the blade dimensions. 

3. Solee equation (36) for the eigenfunctions and determine the 
eigenvalues by application of the boundary conditions on the radial ve- 
locity v. For many applications, F(r) can be taken equal to zero and 
all constants in the series solution for v taken equal to zero except the 
first ones. For such applications, the magnitude of hy will be approxi- 
mately 7. | 

4. Continue the summation procedures that result in relations 
Similar to equations (41) and (42). 


5. Pick u and u,- Notice u, must satisfy the condition that 


r 
oO 


je wee 
r 
i 


where Yr, and r. are the inner and outer boundaries of the cylindrical 
annulus. 

6. Determine w from the first of equations (40). There will 
remain an undetermined constant to choose. 

7. Choose the pressure rise across a stage and the degree of re- 
action desired. The pressure rise across a stator 2P 6 can then be de- 


termined. 
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8. Determine Wo {rom the second of equations (40). The unde- 
termined constant is then evaluated by use of the second of equations (39). 
9. Evaluate the function ¢. 
10. Calculate the axial and whirl velocities for both the stator 
and rotor by using the equations derived in step 4. 


11. Compute the blade angles by use of equations (28) and (29). 


It was decided to use this theory to compare with the blade design 
given in Ref. 1 for solid body blading. The steps will be taken as listed 
above and the results presented here and in Table 1 on page 58. 

l. The loading function was taken to be constant for ease in 
computation. A triangular distribution would be more nearly correct. 
2. The ratio a/b was takentobe 1/3. 
3. F(r) was taken equal to zero and hy found to be approximately 
3.18. The function F(r) was subsequently checked and found negli- 
gible compared to (} fey” : 
4. Equations (41) and (42) were then correct as stated in the 
theory. 
5. wand u » Were then chosen to approximate as closely as pos~ 
sible the axial flow used in the design of the solid body blading in Ref. 1. 
The relations were 
U/U,, = = .560€ + .916 
uo/us = .185€ -.151 
where € is a non-dimensional radius based on the tip radius and Us is 
the rotor tip speed. 


6. Equation (40) was used to determine w in the form 
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w = e - . a e 
(=) 342€-.1578 * 


The constant c, was determined by choosing w=.385 at €=0.8. 
7. The ideal pressure coefficient for the compressor was 0. 40. 


Then by definition 


Td wit ee 9.40 ., 2 
p p 3 i 


where Un, is the blade tip speed that has been used throughout this 


analysis in order to non-dimensionalize the velocity components. Fifty 
per cent reaction in this compressor occurred at € =0.8. Since 


2p, =P, » we obtained the following relation 
2 

ay e 

8. The equation for Wo then became 


Po/p = 0.05 u 





Ww 
Oo _ I 2 . 009 
ar = a (052 6" - 0858 - >). 
T w/e E 


9. The function 6 was evaluated numerically and is listed in 
Table 1. 

10. and ll. The values of u and w are listed in Table 1 for 
different values of €. Similarly, the entering and leaving angles Ye 
Yo? By » and B, are tabulated. The values of these angles used in the 


original design are also tabulated. 





-56- 
VY. CONCLUDING REMARKS 


For most compressor designs occurring in practice, the aspect 
ratio of the blades is neither very low nor very high. It falls some- 
where in between. The theory of Part I of this thesis is applicable for 
blading of low aspect ratio, where it is assumed that equilibrium flow 
patterns exist at the leading and trailing edges of the blade rows. A 
blade row is not affected by neighboring blade rows except through the 
incident velocity distributions. The theory of Part II for blades of high 
aspect ratio contains the influence of disturbances of potential type that 
can be felt upstream as well as downstream. 

The close agreement between the results of the two theories for 
an aspect ratio of three indicates that there the flow between blades is 


very near its equilibrium state at the leading and trailing edges. 
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